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ABSTRACT 

An  analysis  is  made  of  the  distortion  of  a 
plane,  supersonic  free  jet  by  a  strong  magnetic  field. 
The  field  is  created  by  a  straight  wire  placed  outside 
of  the  jet  and  normal  to  the  plane  of  flow.  The  entire 
flow  field  is  obtained  numerically  by  the  method  of 
characteristics  for  several  locations  of  the  wire  when 
the  magnetic  Reynolds  number  is  small.  The  magnetic 
field  tends  to  decelerate  the  flow  to  sonic  velocity 
without  the  formation  of  shocks.  The  resultant  force 
acting  on  the  wire  is  found  to  vary  linearly  with  the 
interaction  parameter  N  for  values  of  N  as  large  as  T^.l. 


I.  INTRODUCTION 


A  partially  ionized  gas  issues  from  an  orifice  with  super¬ 
sonic  velocity  and  atmospheric  pressure.  It  forms  a  two-dimensional 
free  jet  which  would  continue  unimpeded  if  it  did  not  encounter  a 
magnetic  field  since  viscous  and  heat  transfer  effects  are  neglected 
This  report  considers  the  interaction  of  such  an  ionized  jet  with  a 
powerful  magnetic  field. 

Two  cases  are  studied  which  differ  only  in  the  location 
of  the  wire  producing  the  magnetic  field.  In  the  first  case,  shown 
in  Figure  1,  the  wire  is  placed  far  downstream  of  the  jet  exit  so 
that  the  boundary  conditions  of  uniform  flow  at  the  exit  are  still 
valid.  In  the  second  case,  the  line  has  been  moved  to  the  plane  of 
the  exit  of  the  jet  as  shown  in  Figure  2.  Although  the  magnetic 
field  enters  the  nozzle,  the  same  boundary  conditions  are  used  so 
that  comparison  can  be  made  of  the  effect  of  the  location  of  the 
wire.  It  can  be  assumed  for  this  case  that  the  gaB  becomes  conduct¬ 
ing  just  at  the  exit  of  the  nozzle. 

In  order  to  uncouple  Maxwell's  equations  from  the  fluid 

mechanics  equations,  the  magnetic  Reynolds  number  R  is  assumed 

m 

small.  The  interaction  parameter  N  is  assumed  to  be  of  first  order, 
requiring  numerical  solution  of  the  resultant  non-linear  equations 
by  the  method  of  characteristics.  Results  of  a  further  simplifica¬ 
tion  by  linearization  for  small  N  are  presented  in  Reference  1* 
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II.  DERIVATION  OF  THE  STEADY  FLOW  EQUATIONS 

The  steady  flow  of  a  partially  ionized  gas  in  the  presence 
of  a  magnetic  field  is  described  by  the  following  continuum  rela¬ 
tions: 


Momentum 

p(q  »  V)q  +  Vp  s  j  x  B  (1) 

Continuity 

7  .  pq  a  0  (2) 


Energy 

q  •  7  (T  +  =  q  •  (j  x  B)  ♦  (5) 

a 

Ohm's  Law 

j  =  o(E  +  q  x  B)  (4) 


Maxwell ' s  Equations 
7  x  B  =  yj* 

7  x  E  =  0 

7  .  B  =  0 


(5) 

(6) 
(7) 


First  and  Second  Laws  of  Thermodynamics 

T7?  =  Vi-  32  (8) 

P* 

Gas  Law 


P 


pRT 


(9) 
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Vector  quantities  are  indicated  by  an  arrow.  These  in- 
elude  the  velocity  q,  the  current  j,  the  electric  field  E,  and  the 
magnetic  field  B.  Scalar  quantities  in  their  dimensional  form  are 
indicated  by  a  tilda.  These  include  pressure  p",  density  f>,  tempera- 

"w  ew 

ture  T,  enthalpy  i,  entropy  S,  electrical  conductivity  o,  and  per¬ 
meability  p. 

The  inclusion  of  magnetogasdynamic  effects  has  added  new 
equations  as  well  as  additional  terms  to  the  existing  gas  dynamic 
equations.  The  magnetic  body  force  term  has  been  added  to  the  right 
hand  side  of  the  momentum  equation,  and  the  right  hand  side  of  the 
energy  equation  contains  two  additional  terms  which  represent  the 
work  done  by  the  magnetic  body  force  and  Joule  heating.  The  effects 
of  heat  transfer  and  viscous  dissipation  have  been  neglected.  For 
the  two-dimensional  jet  the  electric  field  must  vanish. 

E  =  0  (10) 

Making  use  of  this  fact  to  simplify  Ohm's  Law,  .the  energy  equation 
(3)  reduces  to 

q  •  7  (i  =  0  (11) 

which  states  the  total  enthalpy  H  must  be  conserved.  Thus,  the  work 
done  by  the  magnetic  body  force  just  balances  the  Joule  heating  so 
that  the  flow  remains  isoenergetic .  The  more  general  case  where  E 
does  not  vanish  could  still  be  solved  by  the  method  of  character¬ 
istics,  but  the  choice  of  dependent  variables  would  be  different 
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since  the  total  enthalpy  would  have  to  be  calculated  for  each  point 
on  the  characteristic  net* 

Making  use  of  the  vector  identity 

(q  .  b)q  *  V  -  q  x(V  x  q)  (12) 


the  pressure  can  be  eliminated  between  Eqs.  (l)  and  (8)  leading  to 


VH  -  q  x(V  X  q)  -  ^VS  =  ~  (q  x  B)x  B 

P 


(13) 


Scalar  multiplication  of  this  equation  by  q  results  in 


q  •  VS-  .  i  (J  x  B)2 


(14) 


In  scaler  form,  this  is 


^ 

.  ins 

dx  dy  pT^  y  ^ 


(15) 


From  these  relations  it  is  clear  that  the  flow  will  be  rotational  and 
entropy  will  increase  along  streamlines* 

Using  Eqs*  (8)  and  (9),  the  pressure  can  be  expressed  in 
terms  of  the  other  variables. 


0JL  =  2JSL.  ££  +  ^2  d£ 

v  ^ 

dx  '  dx  dx 

iE.  =  £®_ 

dy*  v  dy  dy 


♦ 


(16) 
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The  final  form  of  the  momentum  Eq.  (l)  ie  obtained  by  elimination  of 
the  pressure  with  Eq.  (16)  and  the  current  with  Eq.  (4). 

+  +  +  =  -  35  fiE  -  vB  }  (17) 

dx  ay  Y  ax  ax  ^  y 

+  +  +  ziQz  . sQ  (18) 

dx  dy  T  dy  dy  *  y 


The  scaler  form  of  the  continuity  equation  is 

a(pu)  +  a(pv)  _  0 

dx  dy 


(19) 


At  this  point  it  is  convenient  to  non-dimensionalize  by 
introducing  the  new  variables 


c 

8 

V 

=  v/u. 

X  X 

T  =  T/T 

oa 

^  / 

*/*» 

P 

*  'P/P. 

B  =B/Bo 

s  •  s/c 

V 

(20) 

ou  LB 

€9  0 

N 

oB2L 

0 

s  1  ■ 

P  u 

r  CD  QD 

The  governing  equations  are  tabulated  on  the  following  page* 
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9u 

dx 

9u 

9y 

9v 

dx 

9v 

9y 

d£ 

dx 

9£ 

dy 

as 

dx 

dS 

ay 

pu 

pv 

0 

0 

2 

a 

0 

2 

£JL_ 

YC 

V 

0  = 

-  NB  (UB  -vB  ) 
y  y  x 

(21) 

0 

0 

pu 

pv 

0 

2 

a 

0 

2 

NB  (UB  -  vB  ) 
x  y  x 

(22) 

P 

0 

0 

P 

u 

V 

0 

0  = 

0 

(23) 

0 

0 

0 

0 

0 

0 

u 

V  e 

-SklJB  -  »B  )2 
p«2  * 

(24) 

dx 

dy 

0 

0 

0 

0 

0 

0  = 

du 

(25) 

0 

0 

dx 

dy 

0 

0 

0 

0  * 

dv 

(26) 

0 

0 

0 

0 

dx 

dy 

0 

0  = 

dp 

(27) 

0 

0 

0 

0 

0 

0 

dx 

dy  « 

dS 

(28) 

In 

addition, 

2 

a  = 

1  ♦ 

X^M2 
2  * 

(1  -  u2 

-v2) 

(29) 

K  ■ 

y(Y-l)Mf 

(30) 

When  the  flow  is  supersonic,  the  set  of  equations  (21)  - 
(28)  is  hyperbolic.  Since  the  set  of  equations  is  linear  in  the  first 
derivatives  of  dependent  variables  u,  v,  p  and  S,  these  derivatives 
can  be  treated  as  unknowns  and  solved  for  in  the  usual  fashion.  Once 
the  first  derivatives  have  been  obtained,  higher  derivatives  can  be 


-  7  - 


calculated  and  used  in  a  Taylor  series  development  near  the  initial 
line  to  g^ve  the  solution  over  a  small  region*  This  could  presumably 
be  repeated  until  the  solution  is  obtained  over  the  entire  region  of 
interest* 

There  are,  however,  exceptional  lines  called  characteristics 
along  which  the  determinant  of  the  coefficients  of  the  left  hand  sides 
of  Eqs.  (21)  -  (28)  vanishes.  In  order  to  have  a  solution,  the  de¬ 
rivatives  must  be  discontinuous  across  the  characteristic.  The 
dependent  variables  must  satisfy  certain  compatibility  relations  ob¬ 
tained  by  replacing  one  of  the  columns  of  the  determinant  by  the  right 
hand  sides  of  Eqs.  (21)  -  (28).  These  relations  are  used  later  to 
obtain  the  values  of  the  dependent  variables  over  the  entire  flow  field. 

The  three  characteristics  for  supersonic  flow  and  their 
corresponding  compatibility  relations  are  given  below. 


Characteristic  I  (streamline) 


dy  v 
dx  ”  u 


(31) 


Compatibility  I-a  (energy) 


.  J2L(uB  -  vB  )2 

dy  pa2  y 


(32) 


Compatibility  I-b  (momentum) 

N(uB  -  vB  )  +  pu  +  pv  4™  +  ^5  4^  ♦  4” 

y  x  K  dx  K  dx  „2  dx  „«2  dx 


M 


ym: 


(33) 


“  8  - 


Characteristic  II  (-  sign),  III  (+  sign)  Mach  lines 


(34) 


Compatibility  II  and  III 


-  N(uB  -  vB  ) 

y  x 


2 

^  .  N  N(y-1)m£[*uB  -  vB  1  ["v  ~  -  ul 

(  g.Bj - J , 

pa 


3  m 

a  T  1  dS  u  du  v  dv~l  u  dy  *  v 
JJ2U(Y-i)  dy  ft2  dy  dy  J  dx  _ 

•  dv 


(35) 


A  typical  characteristic  grid  is  shown  in  Figure  3*  Through  each  of 
the  points  1,  2,  3,  three  characteristics  intersect.  In  the  figure, 
the  streamlines  are  dashed  and  the  Mach  lines  are  solid.  Construction 
of  a  network  of  such  grids  over  the  entire  flow  field  gives  the  dis- 
tribution  of  velocity,  density,  entropy,  and  other  flow  quantities. 

The  details  of  this  are  given  in  a  later  section. 

In  addition  to  the  flow  variables  there  is  another  quantity 
which  is  of  interest.  This  is  the  reaction  force  on  the  wire  which 
produces  the  magnetic  field.  For  each  elemental  area  of  the  jet 
acted  upon  by  the  magnetic  body  force,  there  is  an  equal  and 
opposite  reaction  at  the  wire  according  to  the  Biot-Savart  Law*. 

The  components  of  body  force  per  unit  area  appear  on  the  right  hand 

*  Newton's  Law  does  not  strictly  apply  for  these  magnetic  inter¬ 
actions.  This  point  is  discussed  further  in  Reference  2. 
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sides  of  Eqs.  (21)  and  (22).  Integration  over  the  area  of  the  jet 
followed  by  a  change  In  sign  givee  the  components  of  reaction  foroe 
at  the  wire  as 


(36) 

(37) 


The  magnetic  field  is  produced  by  a  steady  current  flowing 

through  a  straight  wire  located  at  (x  ,  y  ).  The  field  lines  are 

c  c 

concentric  circles  about  the  wire,  and  the  field  strength  varies 
inversely  with  the  distance  from  the  wire.  The  components  of  this 
field,  non-diaensionalised  with  respeot  to  the  field  strength  at 
(xc,  2)  are 


(x  -  x  )  (y  -  2) 
c  c 

(x  -  X  )2  ♦  (y  -  y  )2 
c  c 


(38) 


B 


(yc-y)(yc-2) 

(x-x  )2  ♦  (y -y  )2 
c  c 


(39) 


III.  BOUNDARY  CONDITIONS 

When  the  flow  is  supersonic,  the  system  of  differential 
equations  is  hyperbolic  and  boundary  conditions  must  be  specified 
on  a  U-shaped  boundary.  Initial  conditions  specifying  uniform  flow 
at  the  nozzle  exit  are 

x  =  0:u  =  p  =  l,v  =  S  =  0  (40) 

Along  the  edges  of  the  free  jet  the  pressure  must  be  atmospheric. 

2 

This  means  p  *  1  or  pa  =  1.  Using  Eq.  (29)  for  the  speed  of 
sound , 

p[l  +  ^  Mf  (1  -  u2  -  v2)]  *  1  (41) 

The  shape  of  the  jet  is  not  known  a  priori,  but  must  be  determined  by 
applying  this  boundary  condition  to  the  two  edge  streamlines  as  the 
calculations  proceed  downstream. 
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IV.  NUMERICAL  SOLUTION 

The  three  families  of  characteristics  intersect  at  each 
of  the  infinite  number  of  points  in  the  jet.  In  this  section  a 
numerical  solution  is  obtained  by  using  a  finite  number  of  points  in 
the  jet.  To  start  the  numerical  solution  the  width  of  the  jet  at 
the  nozzle  exit  is  divided  into  a  number  of  equally  spaced  points. 
The  characteristics  through  these  points  are  found,  and  the  con¬ 
struction  of  the  characteristic  mesh  is  continued  downstream  in  the 
usual  stepwise  fashion. 

The  equation  of  the  characteristics  and  their  correspond¬ 
ing  compatibility  relations  are  placed  in  finite  difference  from 
retaining  only  first  order  differences  for  derivatives.  The  other 
quantities  appearing  in  the  equations  are  approximated  by  their 
average  value  along  the  characteristic.  The  finite  difference  form 
of  Eqs.  (31)  -  (35)  for  the  typical  grid  shown  in  Figure  3  is* 

Characteristic  I  (streamline) 

v34(x3  "  x45  ”  u34(y3  "  y4J  *  0  (lf2) 


Compatibility  I-a 


Y(t-DM2N(x  -x4). 
p34a34u34 


v  ”  v34Bx 
y54  *34- 


)2 


(43) 
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Compatibility  I-b 


-*TT  [u34(u3"u4>  +v34(v3"v4)] 
a34 


(44) 


Characteristic  II 


U23V23 


(45) 


Compatibility  Il-a 

[(u23  "  V23)(x3  "  X2)  + 2u23V23(y3  ’  y2>]  U3  +  [2u23V23(x3  ”  X2} 

♦  (v23-u23)(y3-y2)]v3  =  [(U23’V23)(X3*X2) 

+  2u23V23(y3‘y2)]u2+  [2u23(x3-x2)  +(v23’U23)(y3'y2)]V2 


P^U23By23  “  V23Hx23XV23(x3  *  X2}  '  u23(y3  *  y2)IBy23(y3  ‘  y2 


"1  1 
+  B  (x,  -  X-)  -  f'  *-•  ■* —  u„(x,  -  X-)  +  v.,(y,  -  y5) 

x23  5  2  j  r(r-l)Mf  L  23  5  2  23  3  2  J 

Cr-  Dm2n, 


P23a23 


n  up  -ic  -  vc 

riv25<x5  -«2)-  u2J(y5  -  y2>  J  V25°y„  ‘  V25B^j> 


(46) 
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Characteristic  III 


<y3’yl)[(TI)  ““13]  =  (x3  '  X1)[‘U13V13  ^u13  +  v13  ’Off] 


(47) 


Compatibility  Ill-a 


[(u13  '  v13)(x3  '  xl>  *  2u13v13<y3  -  VK  *  [2“i3»i3<*3  * *l> 

*  <’l3'u13><y3"yl,],3  *  [<u13‘v13)<x3‘xl) 

*  2u13v13(y3‘yl)]”l*  [Ju13,13(x3‘xl)  *<,13‘u13)(,3'yl)],l 


*  p7^U13By13  '  V13“x13IV13<x3  -  xl>  -"I3(y>  *yl)IBy,  ,(y3  -yl) 


'13 


1  ~  1 

♦  B  (x,  -  x. )  ^ — 5 —  u,  ,(x,  -  x.  )  ♦  v-  ,(y,  -  y, ) 

X13  5  1  J  r(T-i)M2  L  15  5  1  13  5  1  J 


ivl  3(*  3‘xl)-ul3(y  3'yi)Jvu: 


pl3al3 


2,  n2 
B-  ; 


!3  y15  13  xxy 

(48) 


Point  4  in  Figure  3  lies  on  the  straight  line  connecting 
points  1  and  2  so  that 


yl’y2  y4"y2 


X1  "x2 


x4  ”  x2 


(49) 
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Assuming  thst  the  dependent  variables  vary  linearly  from  1  to  2, 

u4  =  1 


*  ("i-«2>C1-yP 

(50) 

*<’x 

(51) 

*s->2>(^2 

(52) 

(53) 

A  single  subscript  refers  to  a  particular  point  in  Figure  3, 
while  a  double  subscript  indicates  an  average  of  the  values  at  the 
two  points. 

For  the  typical  grid  shown  in  Figure  3,  the  location  of 
points  1  and  2  and  the  values  of  the  dependent  variables  u,  v,  p,  and 
S  at  these  points  is  known.  Using  the  twelve  Eqs.  (42)  -  (53)*  the 
twelve  unknowns  Xy  y y  u y  Vy  Py  S^,  x^,  y^,  u4*  v^,  p4*  and  s4  are 
found.  An  iteration  procedure  is  used  because  u^f  v^,  u^,  and  v^ 
appear  non-linearly.  For  the  first  iteration,  the  average  value  of 
a  quantity  which  involves  point  3  such  as  u^  is  approximated  by  u^. 

The  first  quantities  to  be  computed  are  the  coordinates  Xy 
y^  of  the  intersection  of  the  II  and  III  characteristics.  These  are 
obtained  by  simultaneous  solution  of  Eqs.  (45)  and  (47).  Similarly* 
the  intersection  x^,  y4  of  the  streamline  characteristic  with  the 
straight  line  connecting  points  1  and  2  is  found  from  Eqs.  (42)  and 
(49).  The  values  of  u^,  v4,  p4,  and  S4  can  then  be  found  from 
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Eqs.  (50)  -  (53)«  The  compatibility  equations  (43)  and  (44)  yield 
values  for  and  p^,  respectively,  u^  and  v^  are  obtained  by  simul¬ 
taneous  solution  of  Eqs.  (46)  and  (48).  This  completes  the  first 
iteration*  These  calculations  are  repeated  until  the  numerical 
values  of  the  dependent  variables  converge  to  the  desired  degree  of 
accuracy. 

All  points  in  the  interior  of  the  jet  are  solved  in  the 
manner  described  above.  Points  on  the  upper  and  lower  edges  of  the 
jet  must  be  treated  in  a  different  manner  since  only  two  families 
of  characteristics  intersect.  A  typical  mesh  on  the  lower  boundary 
of  the  jet  is  shown  in  Figure  4.  The  location  of  points  1  and  4  and 
the  values  of  dependent  variables  at  these  points  are  known.  The 
six  unknowns,  Xy  y y  u y  Vy  p ^  and  S^,  are  found  from  Eqs.  (42), 
(43),  (44),  (47),  (48)  and  the  condition 

P3  =  1  (54) 

which  states  that  the  pressure  along  the  edge  streamlines  of  the  jet 
must  be  atmospheric.  A  more  useful  fora  of  this  equation  is 

Pjt1  ♦  ^  -  Uj  -  v*)]  =  l  (55) 

The  calculations  for  a  boundary  point  begin  with  the  simul¬ 
taneous  solution  of  Eqs.  (42)  and  (47)  for  x^  and  y y  Next,  is 
obtained  from  Eq.  (43),  and  p^  from  Eq.  (44).  Then  u^  and  v^  are 
found  by  simultaneous  solution  of  Eqs.  (48)  and  (55).  This  completes 


one  iteration 
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The  iteration  procedure  described  above  did  not  converge 
when  average  values  between  points  1  and  3,  and  points  4  and  3  were 
used  for  u,  v,  and  p.  In  order  to  make  the  iteration  procedure 
converge,  we  introduce  the  weighting  factors  a,  0,  and  £  so  that 


Ui3 

=  *[Ul 

+  auj] 

(56) 

Vi3 

■  «*i 

♦ 

(57) 

pi3 

*  it.i 

+  £p33 

(58) 

The  subscript  i  is  1  or  4  for  the  lower  boundary  of  the  jet  and  is  2 
or  4  for  the  upper  boundary*  These  equations  reduce  to  those  for 
average  values  when  asfisgsl.  For  the  examples  treated  in  the 
report,  rapid  convergence  was  obtained  with  a  «  0.4,  p  =  0*9,  and 
£  =  0.7. 

The  last  quantity  to  be  computed  is  the  reaction  force  at 
the  wire.  This  is  obtained  by  summing  the  components  of  the  equal 
and  opposite  body  force  acting  on  each  diamond-shaped  mesh.  For  the 
typical  mesh  shown  in  Figure  3,  the  area  is 

A  =  /s(s-a)(s-b)  (s-c)  ♦  /s'  (s'-cks’-dMs'-e)  (59) 


s  = 


s '  = 


where 


£(a  +  b  +  c) 
fc(c  +  d  +  e) 


(60) 
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Using  Eqs .  (36)  and  (37),  th«  components  of  the  body  force 


acting  over  the  area  A 

are  approximated  by 

F 

X 

*  »{[uBy  -  »BJ 

(61) 

F 

y 

■  -“{K  -  yBx]  Bxl  .  A 

(62) 

where 

{ 

j  ■*£{). 

avg  i=l  i 

(63) 
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V.  RESULTS  AND  CONCLUSIONS 

The  first  case  to  be  considered  is  shown  in  Figure  6.  The 
jet  leaves  the  nozzle  exit  with  M#  =  2.0  and  interacts  with  the  mag¬ 
netic  field  produced  by  a  wire  located  at  (0,3).  Calculations  were 
made  for  the  interaction  parameter  N  =  0.1.  The  interaction  bends 
the  jet  upward  and  decelerates  the  flow.  At  X  =  10,  the  upper  edge 
of  the  jet  has  been  displaced  twice  as  much  as  the  lower  edge  so 
that  the  width  of  the  jet  has  doubled.  Contour  lines  of  constant 
Mach  number  and  u  appear  in  Figures  7  and  8,  respectively.  Since 
the  flow  is  rotational,  the  entropy  increases  along  streamlines  as 
indicated  by  Figure  9*  The  density  decreases  in  the  streamwise 
direction  even  though  the  velocity  is  also  decreasing.  Contour  lines 
of  constant  density  are  plotted  in  Figure  10. 

In  Figure  11  the  components  F  and  -  F  of  the  reaction 

^  y 

force  at  the  wire  are  plotted  for  two  locations  of  the  wire,  the 

solid  curves  for  (0,3)  and  the  dashed  curves  for  (0,1.5) •  For  the 

range  of  N  considered  the  components  vary  linearly  with  N  except 

near  N  =  0.  The  magnitude  of  the  reaction  force  increases  directly 

with  N  and  its  inclination  is  fixed  in  the  downstream  direction 

toward  the  jet.  Thus,  the  jet  tends  to  drag  the  wire  with  it.  Care 

must  be  taken  in  comparing  the  two  sets  of  curves  for  the  same  value 

of  N,  for  N  is  proportional  to  Bq,  the  value  of  the  field  at  the 

origin.  To  keep  Bq  fixed,  the  wire  at  y^  =  3  must  carry  twice  the 

current  of  the  wire  at  y  =  1»5»  Although  at  y  =  1.5  the  wire  is 

c  c 

closer  to  the  jet  near  the  exit,  the  downstream  portion  of  the  jet 


-  19 


is  about  equidistant  from  both  locations  due  to  the  bending  of  the 

jet*  These  factors,  along  with  the  greater  area  downstream  for 

interaction,  leads  to  the  greater  force  on  the  wire  at  the  outer 

location  y  =  3*0. 
c 

The  next  two  figures  show  the  jet  when  the  wire  is  located 
at  (20,3)  which  is  far  downstream  of  the  nozzle  exit  at  x  =  0.  In 
Figure  12  M  =  4.0  and  N  =  0.1.  The  jet  slows  down  and  reaches  the 
sonic  velocity  near  x  =  16.  The  two  edge  streamlines  are  almost 
symmetric  and  bend  so  that  the  width  of  the  jet  at  x  =  16  is  2.5 
times  the  exit  width.  It  is  interesting  to  note  the  smooth  decelera¬ 
tion  without  the  appearance  of  shocks  as  in  the  case  of  the  channel 
flow  considered  in  Reference  3* 

In  Figure  13  Mm  =  2.0  and  N  =  0.05*  The  interaction  is 
weak  enough  to  permit  the  jet  to  pass  the  wire  at  x  *  20.  Along  the 
upper  streamline  at  x  »  27  the  Mach  number  has  reached  unity.  The 
upper  streamline  has  been  displaced  very  little,  while  the  lower 
streamline  has  moved  downward  to  increase  the  width  of  the  jet  by  a 
factor  of  two.  The  components  of  body  force  are  F^  =  0.22  and 
F  =  -  0.012. 
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Figure  1.  INTERACTION  OF  A  PLANE  JET  WITH  A  MAGNETIC  FIELD.  WIRE 
LOCATED  DOWNSTREAM  OF  THE  JET  EXIT. 
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EPICAL  CHARACTERISTIC 


0RID  F0R  interior  or  JET. 


Figure  4.  TYPICAL  CHARACTERISTIC  GRID  ALONG  LOWER  BOUNDARY  OF  JET 
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Figure  8.  CONTOUR  LINES  OF  CONSTANT  u. 


Figure  9.  CONTOUR  LINES  OF  CONSTANT  s 
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E  AT  THE  WIRE 


Figure  12.  JET  SHAPE  FOR  M  =  4.0  AND  N  =  0.1 


Figure  13.  JET  SHAPE  FOR  Ma  =  2.0  AND  N  =  0.05 


